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Abstract 

^\ Complex networks are all around us, and they can be generated by 

simple mechanisms. Understanding what kinds of networks can be pro- 
duced by following simple rules is therefore of great importance. We 
investigate this issue by studying the dynamics of extremely simple sys- 
terns where are 'writer' moves around a network, and modifies it in a way 
that depends upon the writer's surroundings. Each vertex in the network 
has three edges incident upon it, which are colored red, blue and green. 
This edge coloring is done to provide a way for the writer to orient its 
!_h movement. We explore the dynamics of a space of 3888 of these 'colored 

trinet automata' systems. We find a large variety of behaviour, ranging 
from the very simple to the very complex. We also discover simple rules 
that generate forms which are remarkably similar to a wide range of nat- 
ural objects. We study our systems using simulations (with appropriate 
visualization techniques) and analyze selected rules mathematically. We 
arrive at an empirical classification scheme which reveals a lot about the 
kinds of dynamics and networks that can be generated by these systems. 
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1 Introduction 



Many advances in complex systems research have come from identifying sim- 
ple deterministic systems which can generate complex behavior. For example, 
cellular automata [l] are used to model many systems in physics and biology, 
string rewrite [2] systems are used to model plant development, and chaotic 
differential equations [3] are used in many areas. The most popular models 
of complex network growth are probabilistic 4,5. These models have yielded 
many insights, but they do not offer an explanation as to where the complexity 
in networks comes from. The complexity in probabilistic models comes from 
the randomizing mechanisms they arc based on. 

In this paper we study a class of simple deterministic network growth models. 
Our systems can produce a wide variety of behaviour, ranging from the very 
simple to the very complex. They offer insights into how complex structures 
can be generated, because each aspect of their behaviour can be traced back 
to a deterministic cause. The simple nature of our models allows us to quickly 
generate exact pictures of their dynamics. This gives us the ability to easily 
explore the behaviour of large numbers of systems. By taking an unbiased look 
at the dynamics of many simple rules we see what types of network growth are 
easily generated by simple computational processes. 

In Turing machines [7] a writer moves along a one dimensional tape, and 
rewrites symbols using local rules. Our network automata systems are like 
a generalization of this idea -where the writer moves around a network, and 
rewrites it on a local level. The way the writer moves and modifies its local 
structure is determined by its surroundings. The network the writer runs on is 
a trinet (i.e., each vertex has three connections). Trinets are also known as cubic 
graphs and three-regular graphs. There are many natural examples of trinets, 
such as two dimensional foams and polygonal networks formed by biological cells 
and cracks in the earth. We focus on trinets because they are easy to manipulate 
and expressive. Any network can be represented as a trinet by replacing vertices 
that have more than three connections with "roundabout-like" circles of vertices 
with three connections [191. 



One of our goals is to find simple deterministic rules that generate complex 
dynamics. An obstacle to this goal is the inability of deterministic rules to break 
symmetries. For example, if the writer is at one of the corners of a cube then 
it has no deterministic way to select which of its "similar looking" neighbors it 
should move to next. We circumvent this obstacle by supposing the edges of 
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our trinet are colored red, blue and green in such a way that touching edges 
have different colors (such an edge coloring is also known as a Tait coloring [6]). 
This allows us to specify how the writer should move by stating which colored 
edge it should move along in different situations^ 

Our main goal is to understand about what kinds of behaviour can be gener- 
ated by colored trinet automata. To achieve this, we do a thorough exploration 
of the dynamics of a space of simple rules. Using pictures and analytic tech- 
niques we sort these rules into three classes (fixed points, repetitive growth and 
elaborate growth) according to their behaviour. We discuss these classes in 
sections |2|3| and [4] In section [5] we exhibit more general rules, which can pro- 
duce other exotic behaviour such as persistent complex behaviour, periodically 
changing networks and production of hyperbolic space. In section [5] we also 
present many rules that can produce networks that look remarkably similar to 
real world objects. 

We have tried to keep the main text as free from equations as possible. This 
was done to increase readership, and because complex phenomena are often 
better explained using pictures. We also studied our systems mathematically, 
and the appendix contains proofs which will give the mathematically inclined 
reader significant insights into some of our system's dynamics. 



1.1 Related Literature 

Many interesting deterministic network growth models have been previously 
considered. These include the growth models discussed in [8] and the biologically 
inspired models considered in [9 -12 , which can produce rich and complicated 



self replicating structures, despite extremely simple rules. Self replication in 



adaptive network models was also considered in 13 . In each of these cases, 
every part of the network gets updated in parallel. In our systems, by contrast, 
only a small piece of the network gets updated at any time. This makes our 
systems easier to implement because there is no need for distant parts of the 
network to have synchronized clocks. It is true that many real world networks 
grow in parallel, however evolving our systems for many time steps often leads to 
structural modifications that are equivalent to global/parrelel rewrite operations 



(see subsections 4.1 and 4.4 1. 



1 We could remove the reliance of edge color in our systems, by replacing the red, blue and 
green edges with different uncolored structures, and modifying our rewrite rules accordingly. 
Although this would make the systems look more complicated. 
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Our work was inspired by the work of Tommaso Bolognesi 19 , who studies 
the dynamics of planar trinet automata. Like our models, these systems involve 
a writer which moves around trinet, making structural modifications as it goes. 
Unlike our systems, the models considered in 19 circumvents the symmetry 
breaking problem by supposing the trinets are embedded in the plane. The 
way the writer behaves in these systems is governed by the structure of the 
faces formed by the planar embedding. Although these models arc fascinating, 
the fact that the dynamics depend on how the network is embedded makes it 
difficult to get a full picture of what is going on at any given time. Our approach 
(of considered edge-colored trinets, rather than planar embedded trinets) allows 
us to easily picture the complete dynamics (and rules) behind the systems we 
consider. 



The experimental approach we use was pioneered by Stephen Wolfram 14 



who uses simulations to reveal a vast array of simple programs that can generate 
complex dynamics. Wolfram finds many network growth models that can gen- 
erate complex dynamics and explores the idea that the physical universe could 
be generated by a simple network growth model. The idea that the universe can 
be generated in exact detail, by a simple program, is one of the most interesting 
conjectures from digital physics 15 16 18 . Wolfram suggests that the correct 
model could be found by doing an automated search of the simplest possibil- 
ities 17 . In order to do such universe hunting one needs intuition about the 



kinds of things that simple adaptive network mechanisms can do. One of our 
aims is to provide some of this intuition. 



1.2 How Our Systems Work 

A colored trinet automata is a dynamical system where a writer moves around 
the vertices of an edge-colored trinet; applying rewrite rules as it goes. Every 
time step some rewrite operation is applied about the writer's current location 
and then the writer moves. We start by considering only extremely simple 
rewrite operations, where the writer's current vertex is either replaced with a 
triangle or left unaltered. The action the writer takes on a given time step is 
determined by the colors of the edges interlinking its neighbours. 

The rules behind a colored trinet automata specify how the writer should 
move, and modify the network in response to its surroundings (i.e., the colors of 
the edges interlinking the writer's neighbors). We picture the rules at the top 
of our figures, and the writer is represented by a black vertex. For example, the 
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Figure 1: An example colored trinet automata. At the top we show the rules of 
the system by indicating how the writer (the black vertex) moves and modifies 
the network in response to its surroundings. Underneath we show the system 
evolving for four updates, starting from the cube. The system reaches a fixed 
point after three updates. 



rule for the system shown in Figure [T] can be described as follows; 

1. If there arc no edges linking the writer's neighbours then the writer then 
moves along a blue edge, and the writer's previous location is replaced 
with a triangle. 

2. If there is exactly one edge linking the writer's neighbours, which is red, 
then take no action. 

3. If there is exactly one edge linking the writer's neighbours, which is blue, 
then the writer along a red edge and the network is left unaltered. 

4. If there is exactly one edge linking the writer's neighbours, which is green, 
then take no action. 

Notice how these instructions correspond to the pictures in the four boxes at 
the top of Figure [T] In general, a rule is just a specification of which struc- 
tural modifications and movements the writer should perform in response to 
the four types of surroundings (depending on colors of edges interlinking the 
writer's neighbors) that the writer can have when there is no more than one 
edge linking its neighbors. Our rules do not specify which actions the writer 
should take in other situations, where there are two or more edges linking its 
neighbors. In these cases we suppose that no action is taken. Actually, the rules 
we consider never generate vertices with two or more linked neighbours and so 
these situations never occur anyway. 
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Figure 2: An illustration of the space of rules we explored. A rule is specified 
by associating each of the four different types of surroundings with one of the 
images to its right. 

1.3 The Rule Space 

In each of our investigations we used the cube (shown at the bottom left of 
Figure [I]) as our initial condition. We studied the dynamics of the space of 
12 x 18 x 18 = 3888 rules depicted in Figure [2] 

These are all the rules which involve possible triangle-replacement, and 
movement along two or less links, with two exceptions. Our first exception 
is that we always replace a vertex with no interlinked neighbours with a trian- 
gle. We make this exception so that each of our rules is guaranteed to change 
our initial cube network. Our second exception is that we never take action 
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when there is a green edge linking the writer's neighbors. We make this excep- 
tion to reduce the number of rules we must consider. If we drop this second 
constraint we find new kinds of behaviour, but the rule space becomes too large 
to thoroughly explore. 

By severely limiting the kinds of operations our trinet automata can perform, 
we reduce the rule space to a manageable size and insure that it contains minimal 
systems that generate certain types of behaviour. We sort our 3888 rules into 
three classes (fixed points, repetitive growth and elaborate growth) according 
to the long term dynamics they generate, starting from the cube. The following 
three sections are devoted to describing the behaviour of the rules in these three 
classes. 

2 Fixed Points 

We say a system has reached a fixed point, where there comes a time after 
which the network no longer changes. 2918 (about 75 percent) of the our rules 
eventually reach a fixed point. Our example system from Figure [I] reaches a 
fixed point after three updates. One reason such a large number of rules end up 
at a fixed point is that our rule space is such that whenever a green edge links a 
pair of the writer's neighbors a system becomes fixed. Indeed 2562 of our rules 
halt their evolution because of this effect. These include the rule which grows 
the largest fixed structure, shown in Figure [3j 

In addition to the fixed points where the writer halts, there are dynamic- 
writer fixed point where the writer continues to move forever, even after the 
network has become static. 162 of the rules in class 1 evolve into dynamic- 
writer fixed points. In 130 of these rules the writer ends up doing a period two 
orbit (see Figure [4]), in the others the writer ends up doing a period four orbit. 

3 Repetitive Growth 

The next most complex type of behaviour observed is repetitive growth. This 
is characterized by the feature that the structure continues to grow, while the 
writer is "trapped" within a particular region of the network -with the form 
of its surroundings changing periodically. Repetitive growth occurs because 
writer's surroundings induce it to generate more structure around itself of the 
same form. 840 of our rules eventually generate repetitive growth (see Figure 
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Figure 4: Of all our rules, this one takes the longest (34 updates) to reach a 
static network (which has 52 vertices). The writer continues to move in the 
period two orbit shown after this network has been generated. 
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Figure 5: A rule that generates simple repetitive growth evolving over the first 
four updates. Once the writer's surroundings are such that a pair of its neighbors 
are linked by a blue edge, the writer performs an update which leads to similar 
surroundings occurring again on the next time step. 

Let us define repetitive growth more precisely. Period p repetitive growth 
is occurring at time t when the network grows arbitrarily large eventually, and 
there is a distance r such that the network on vertices within a distance r of 
the writers position at time t looks identical to the network on vertices within 
a distance r of the writer at time t + p, and the writer never moves to an old 
vertex more than a distance r — 1 from where it was at time step t during the 
interval [t,t + p\- The rule shown in Figure [5] falls into period one repetitive 
growth because after the first update the structure within a distance one of the 
writer always looks similar. Evolving this rule for t > time steps results in a 
network with 8 + 2t vertices and a single triangle on the end of a long "ladder 
like" substructure. 

Although one may test whether dynamics fit our definition using a computer, 
there are more practical ways to spot repetitive growth (see Figure [6]). Networks 
undergoing repetitive growth tend to have an elongated linear, or circular shape 
because they are composed of a series of repeating substructures. Plotting the 
index of the writer over time is an excellent way to see dynamics. However it 
does depend on the way the vertices are indexed within the computer program, 
and this inevitably depends on more than just the pure topological dynamics of 
the system. In our case, when a vertex with index v in an L vertex network is 
replaced with a triangle, we give the vertex of this triangle with a red external 
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Figure 6: On the left we show the networks generated on time steps 6, 7 and 
8 (reading downwards), when evolving from a cube. On the right we plot the 
index of the writer's position over time. Notice how the differences between 
the indices of the writer's position forms a periodic sequence. This is because, 
although network keeps growing, the writer keeps moving in the same relative 
way. This system has period two repetitive growth because the network induced 
on vertices within a distance two of the writer on time step 6 looks identical to 
the induced on vertices within a distance two of the writer at time step 6 + 2. 



edge an index v, and we give the other two vertices the other vertices of this 
triangle, with green and blue external edges, indices L+X and L+2, respectively. 



3.1 Repetitive growth with long transients 

All but four of our repetitive growth rules settle into low (i.e., less than 12) 
period repetitive growth quickly (in less than 40 time steps). The other four 
rules eventually produce repetitive growth (see Figure [7]), although it would 
perhaps be best to describe their behaviour as complex because the systems 
have extremely long transients within which the structures appear to grow 
in a pseudo-random way (this is reminiscent of elementary cellular automata 
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Figure 7: This system has complex behaviour for the first 4994 time steps, and 
then settles into period 454 repetitive growth. We show the rule at the top. On 
the left we show an uncolored plot of the network present on time step 7000. 
The large "handle" developing on the right of this picture is a symptom of the 
repetitive growth. The plot on the right shows the writer index over the first 
7000 time steps. 



11 



110 [14]). Interestingly, evolving the rule shown in Figure [7] from other small 
networks produces different behaviour. The network known as K 3 3 is obtained 
by taking two clusters of three vertices, and linking each vertex in one cluster 
to each vertex in the other. Initiating this system from (instead of the 

cube) leads to dynamics which (again) eventually settle in period 454 repetitive 
growth, although this time the transient lasts for 29964 time steps. Another of 
the four rules which induce repetitive growth after a long transient is symmet- 
rically equivalent to the rule shown in Figure [7J because it can be transformed 
into it by swapping the roles of the red and blue colors. 

The other two rules which induce repetitive growth, after a long transient, 
are not symmetrically equivalent to the system shown in Figure [7J although 
they are symmetrically equivalent to each other. One of these rules is shown in 
Figure [HJ 

The ability of these systems to alter their behaviour after large amounts of 
time leads to the production of heterogeneous substructures. The fact that these 
rules produce a complex "ball" of network followed by a long one dimensional 
structure is reminiscent of the way plants grow by complicating the structure 
around the initial seed and then growing out a long stalk (see the figures from 
section [5]) . 

4 Elaborate Growth 

The remaining 130 rules produce elaborate patterns of growth, leading to self 
similar networks. It turns out that in each of these cases, the way the writer 
moves is effectively one dimensional. The writer repeatedly traverses a one 
dimensional "track" , applying rewrite rules as it goes, and lengthening the track 
with each traversal. In most cases, the fact that the writer is confined to a one 
dimensional substructure can be inferred directly from the rules. In the other 
cases this behaviour is revealed by using appropriate visualization techniques. 

The writer can move backwards and forwards on the one dimensional track 
it is confined to. Although the writer can move in complex ways we can sort the 
rules into two subclasses according to the general nature of the writer's move- 
ment. Either the writer moves around and around a closed loop, in which case 
we say the rule has cyclic writer movement, or the writer moves backwards and 
forwards over line, in which case we say the rule has bouncing writer movement. 
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Figure 8: On the left we show the network obtained after evolving this rule for 
300 time steps. The black arrows track the course that the writer has taken 
over the last 150 time steps (note that the positioning of these arrows depends 
on how we indexed our vertices). The writer movement looks unpredictable 
during this early transience. On the right we plot the writer's index over time. 
The transient behaviour lasts 17615 time steps, then the system settles into 
period 1355 repetitive growth. Similar behaviour ensues when evolving from 
other initial conditions. 
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Figure 9: One of the simplest rules with cyclic writer movement, evolving for 
seven updates. 



4.1 A simple case with cyclic writer movement 

104 of the rules with elaborate growth have cyclic writer movement. We show 
one of the simplest in Figure [9| By looking at the rules of this system one can 
see that the writer must move in a one dimensional fashion. Whenever there 
are no red or green edges interlinking the writer's neighbors (i.e., whenever the 
system is not at a fixed point), the writer moves along a red edge and then a 
blue edge, and the writer's previous location is replaced with a triangle. The 
fact that the writer's movement can be expressed as a combination of steps 
along edges of only two different colors implies that the writer's movement is 
effectively one dimensional. In particular, since any sufficiently long path along 
edges of alternating red/blue color must eventually return to its starting point 
the writer must always remain confined to a cyclic track. In this case the cyclic 



track starts out as the inner face of our cube (see Figure 10 ). On each update the 
writer replaces its current vertex with a triangle and moves two edges clockwise 
around the developing track of edges with alternating red/blue colors. Theorem 
[T] describes the global rewrite operation which the writer effectively performs 
with each complete traversal of its cyclic track. The amount of time successive 
traversals take doubles. 
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Figure 10: An alternative "one dimensional" representation of the dynamics 
of the rule shown in Figure [9] Successive rows, reading downwards, show the 
relevant part of the network on successive time steps. Here the relevant part of 
the network consists of network induced on the vertices which can be reached 
by moving along red or blue edges from the writer. We do not show green 
edges which are not part of triangles since these have no effect on dynamics. 
Each time step the writer replaces its current vertex with a triangle (effectively 
adding a red edge and a blue edge below a green arc) and moves two edges to 
the right. 
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Figure 11: The network generated by the system with rules shown in Figure [9] 
on the first few time steps of the form 2™ — 2 : n > 2. Between times of this form 
the writer makes a complete traversal of its cyclic path and effectively performs 
the rewrite operation shown at the bottom, globally. 



Theorem 1 The way the rule shown in Figure [$| evolves (with the cube as 
the initial condition) is such that for each n > 2, the network present on the 
(2™ +1 - 2)th time step can be obtained by taking the network present on the 
(2™ — 2)th time step and then simultaneously replacing each vertex, with a red 
or blue edge interconnecting its neighbors, with a triangle (see Figure 11). 

Proofs are in the appendix. 



4.2 A rule related to the golden ratio 



In Figure [12] we show a rule with cyclic writer movement with a complicated 
looking growth rate. It turns out that the growth rate can be described exactly 
in terms of the golden ratio 
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Theorem 2 The number of vertices in the network obtained by evolving the 



rule shown in Figure 12 for t > time steps (starting from the cube) is 
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It is pleasing that this seemingly complex growth rate can be described 
simply in terms of the golden ratio, because this number appears in so many 
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Figure 12: On the left we show the network obtained by evolving this rule for 
62 time steps. The number of vertices in the network at time step t grows 
approximately linearly with t. However the way the number of vertices deviates 
from its best linear fit looks complicated (as shown on the right). 
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interesting places in the natural world. This colored trinet automata has similar 
qualitative behaviour to the one considered in subsection |4.1| Once again the 
writer goes around and around a cyclic track consisting of edges with alternating 
red/blue color. Again the system can be reduced to a one dimensional string 
rewrite system. The proof to Theorem [2] is based on relating this system to the 
binary rewrite system with rules — > 01, 1 — > Oil. 



4.3 A rule with complex behaviour 

The rule with cyclic writer movement shown in Figure [13] has a very compli- 
cated looking growth rate. Although this system can be transformed into a 
one dimensional rewrite system (in a similar way to the systems considered in 
subsections 4.1 and 4.2 ) we have not been able to derive a formula for its growth 
rate. This means this system, together with its equivalent red/blue reflection; 
stand out as the most complex rules in our entire seiQ We are unable to pre- 
dict the long term dynamics of these systems (although it appears the pseudo 
random growth pattern continues forever). 



4.4 A simple case with bouncing writer movement 

The 26 remaining rules with elaborate growth exhibit bouncing writer move- 
ment. In these rules the writer moves forwards and backwards along a one 
dimensional track, which grows with successive traversals. The rule in this class 
with the simplest behaviour is shown in Figure [14] This rule can be represented 



in a one dimensional manner similar to Figure 10 except that this time the 
writer is confined to moving on the red and green edges. 

It appears that the writer effectively performs a global rewrite operation 
every time it makes a traversal of the one dimensional track it is confined to. 
Simulations suggest the that the network present on time step 2 n+3 + n (where 
n > 0) can be obtained by taking the network present at time step 2™ +2 + n — 1 
and then simultaneously replacing each vertex with a triangle that has no red 
or green edges interlinking its neighbors and is not part of the external face (see 
Figure 15). 

Not all of the rules with elaborate growth can be reduced to one dimensional 
rewrite systems as directly as the examples we have considered, although the 



2 All of the rules which generate fixed points and repetitive growth have trivial long term 
behaviour. Plotting the writer index over time reveals significant regularities in all the other 
rules exhibiting elaborate growth 
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Figure 13: On the left we show the structure obtained by evolving this rule 
for 100 time steps. The black arrows track the positions of the writer over 
the previous 45 time steps. On the right we show how the growth rate of the 
number of vertices deviates from its best linear fit. This plot reveals considerable 
complexity. 
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Figure 14: A rule with bouncing writer movement evolving for eight updates. 
The writer travels around the red and green edges of the track that starts out 
as the central face of the cube. The writer keeps bouncing off the red edge at 
the bottom of this face and reversing its direction of movement. 
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Figure 15: The network generated by the system with rules shown in Figure 14 
on the first few time steps of the form 2 n+2 — n — 1 : n > 0. Between times 
of this form the writer makes a complete traversal of its linear path and both 
rewrite operations shown at the bottom, globally. 



effectively one dimensional nature of the writer's movement can be revealed by 
plotting the trail followed by writer over several time steps (as in the left of 
Figure 13). 



5 More general rules 

The behaviour of colored trinets with more general rules does not always fall 
into the classes listed above. The set of 3888 rules we enumerated, and discussed 
above, did not allow the writer to take any action when the edge between its 
neighbors is green. If we remove this restriction we can find rules with four active 
parts such as the one shown in Figure |16| In this rule the writer continues to 
move in a random looking way (that is not "one dimensional" as in the rules 
with elaborate growth discussed above) for at least the first 100000 time steps. 
It is an open question whether this rule eventually settles into repetitive growth. 

We can also consider rules which include more general kinds of rewrite op- 
erations, such as the one shown in Figure |17| In this system triangles can be 
replaced with vertices. Unlike all of the systems in our previously considered 
rule set (which either has approximately linear growth rates or reached fixed 
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Figure 16: A rule with four active parts that has complex behaviour. On the 
left we show the structure obtained by evolving this rule for 100 time steps. The 
black arrows track the positions of the writer over the previous 45 time steps. 
On the right we plot the index of the writer over the first 1000 time steps. 
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Figure 17: A system with more general rules that may involve replacing a 
triangle with a vertex. On the left we show the network present on time step 
1000. On the right we plot the number of vertices over time. 



points) the number of vertices in this case grows sub-linearly over time. 

When exploring these more general kinds of rules, one finds many cases with 
sub-linear growth rates that exhibit repetitive or elaborate growth patterns. 
In many cases with sub-linear growth the number of vertices grows like the 
square root of the number of time steps elapsed. Other rules involving triangle 
shrinkage can lead to networks with shapes that oscillate over time. A trivial 
example is shown in Figure 18 (although rules exist with much higher period 
oscillations) . 

We can also consider rules including the so called exchange operation, which 
rewires an edge 13 . In Figure 19 we show a rule which uses this type of ex- 



change operation. The rule depicted in Figure 19 has some interesting spatial 
properties. For a vertex v, the number of vertices a distance k from v grows 
exponentially with k. This suggests the rule generates hyperbolic (i.e., nega- 
tively curved) space. Indeed, this can be verified using the techniques described 
in [22] and |23 . According to these works, a connected network G with vertex 
set V(G) is scaled Gromov hyperbolic (which means it is negatively curved on a 
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Figure 18: In this system the network shape changes periodically. 





Figure 19: A system which includes the exchange rule, which effectively replaces 
a triangle with a square. We show the rule at the top. The structure is generated 
by running this system for 500 time steps. 
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large scale) when gf^jy < §• Here 

6(G) — max{min{d(x, y)+d(y, z)+d(z, x) : x € [b, c], y e [a,c], z G [a, 6]}, a, 6, c € V(G)}. 

Where [i, j] denotes the shortest path from i to j , and c?(i, j) denotes the distance 
from i to j. Also dia(G) = m&x{d(u,v) : u,v £ V(G)} is the diameter of G. In 
our case, the network G obtained by evolving the rule shown in Figure[T9]for 100 
time steps has gf^k = 1.28571, and so it is indeed scaled Gromov hyperbolic. 

The 3888 rules in the set that we initially focused on do not generate neg- 
atively curved space. This can be seen by noting the each of these rules gen- 
erates planar networks, and planar networks have non-negative combinatorial 
curvature]^] (see the Gauss-Bonnet formula 24 ) . These rules preserve planarity 



because the only structural modification they allow is the replacement of a 
vertex with a triangle (and this can be thought of as chopping the corner off 
a polyhedron -a planarity preserving operation). On the other hand, the ex- 
change operations within the rule depicted in Figure |19| allow it to generate a 
non-planar network, when initiated from the cube. 

5.1 The resemblance with real world objects 

It is remarkable how varied the forms produced by trinet automata with very 
similar rules can be. What is even more remarkable is that pictures similar to a 
diverse range of life forms and objects can easily be produced. In Figures [2T| [22} 
23(24] and [25 we show networks (produced by evolving colored trinet automata, 



from the cube) next to real world objects. In many cases there is significant 
similarity. When looking at these pictures it is important to keep in mind that 
a trinet may be represented by a wide variety of pictures, because the only 
relevant information is the connectivity. In other words, one could draw the 
colored trinets in a different way, and their resemblance to these objects would 
disappear. Never the less, each of these trinet pictures where produced using 
Mathematica's standard, spring embedding or spring electrical embedding graph 
drawing algorithm. This means the way the networks are drawn was chosen by 
the plotting algorithm to reflect the topology of the network. In this regard, 
one can argue that there is similarity between the real objects shown and the 
colored trinets pictured because the most natural ways to draw the networks 



3 Combinatorial curvature is a local measure, but the fact none of these networks have 
negative curvature on the large scale was verified by checking that each network which grows 
arbitrarily large has a well defined node shell dimension [16] less than 3. 
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h-hX-Xk-Xl-l 




Figure 20: On the left we show the networks generated by evolving this system 
for 150, 200 and 250 time steps, evolving from the cube. On the right we show 
a photo depicting the development stages of a mushroom, taken from 25 . The 



cap is the first part of the mushroom to develop, and then a long stem grows 
our from the spore at the bottom. Our system generates a ball like network 
initially, and then grows out a linear structure. These pictures could be made 
to look more alike by manipulating the network embedding. The networks are 
drawn using Mathematica's default Spring Electrical Embedding 
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Figure 21: On the left we show the network obtained from evolving our rule for 
46 time steps, starting from the cube. The network drawn using mathematica's 
standard spring embedding algorithm. On the right we show a photo of a turtle 



shell taken from 26 



yield pictures that look like the objects. Figures 21 22 23|24 and 25 also show 
topological similarities between the colored trinets and the real objects. 



6 Conclusion 

We have explored many aspects of the behavior of colored trinet automata. 
We found that the systems in our initial rule space could be sorted into three 
classes - fixed, repetitive growth and elaborate growth. We have shown that 
each class can be described in terms of writer movement. When we look at the 
dynamics of more general rules we find other types of behavior. In particular, 
we find rules which generate sub-linear growth, persistent complex behaviour, 
periodically changing networks and hyperbolic space. We have also noted that 
many trinets with more general rules can produce forms remarkably similar to 
a diverse range of real world objects. 

As minimal models capable of growing complex structure, these systems 
should have some applications. The systems could be realized by creating a 
robot which pulls itself along chains, ropes or silk and drops lines (like a spider) 
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Figure 22: On the left we show the network obtained from evolving our rule 
for 34 time steps, starting from the cube. The network drawn using mathemat- 
ica's standard spring embedding algorithm. On the right we show a photo of a 



caterpillar taken from 27 



to triangulate vertices. This type of realization could perhaps be useful for 
weaving or construction work. The ability of these systems to produce exotic 
network structures with relative case could make them useful in network design. 
In existential graph theory there are many open questions as to whether large 
trinets exist with particular properties and it will be interesting to see whether 
our systems can resolve any of these issues. 

It will also be interesting to see if these systems have any applications in 
modeling. The way structures grow under repetitive growth is reminiscent of 



the way plants grow (see Figure 20 1 . Many of the networks produced by rules 
with elaborate growth rules seem reminiscent of polygonal networks in cracked 
earth or foams. The general idea of having a single writer that rewires a complex 
network is reminiscent of brains and databases where the connectivity is altered 
to store information. 

There many directions that this work can be taken in the future. We have 
explored the dynamics of many colored trinet automata and identified many 
kinds of behavior, however a more general classification scheme is evidently 
required to deal with the systems we described in section [5j Exploring more 
general rules will be exciting, and will surely yield systems with other interesting 
kinds of behavior and applications. There are also many unanswered questions 
about the systems presented here, for example: is there a simple formula for the 



growth rate of system shown in Figure 13 ? Does the complex behavior in Figure 
16 persist forever ? How can rules be characterized according to the structural 




Figure 23: At the top we show the rule. In the middle we show the network 
obtained from evolving our rule for 1600 time steps, starting from the cube. 
The network drawn using mathematica's standard spring electrical embedding 
algorithm (we have not shown edge colors) . At the bottom we show a photo of 
sliced Telescopium seashells taken from [28], 
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Figure 24: On the left we show the network obtained from evolving our rule 
for 1430 time steps, starting from the cube. The network drawn using mathc- 
matica's standard spring embedding algorithm (we have not shown edge colors). 
On the right we show a photo of a coronal slice through a human brain taken 



from 29 











I-A. 









Figure 25: On the left we show the network obtained from evolving our rule for 
100 time steps, starting from the cube. The network drawn using mathematica's 
standard spring electrical embedding algorithm. On the right we show a photo 



of a Necklace made from crochet lace, pearls, and sterling silver taken from 30 
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properties of the networks they produce ? 
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A Appendix 

A word w with alphabet A is a sequence wo...w m -i of m > elements Wi € A 
from w. The length \w\ of such a word is the number elements in it (m in this 
case). The empty word is a word of length 0. Notice that we index the elements 
in our words w from to \w\ — 1. Let A* denote the set of all words with 
alphabet A. Let 

WX = Wo...W\ w \_ 1 X ...X\ x \_ 1 (1) 

denote the word formed by concatenating words w and x. If k > then w 
denotes the word obtained by concatenating k copies of w together (e.g. w 3 = 
www). Also w° is the empty word. 

We let r, b and g denote the colors red, blue and green, respectively. We 
define the type p(v)c of a vertex v in a colored trinet G such that, if none of 
w's neighbours are interlinked then p(v) G — and if exactly one pair of w's 
neighbors are interlinked, by an edge of color c e {r, b, g} then p(v)c = c. 

We let (X*,G*) denote the state of a system at time step t (i.e., after t 
updates, staring from the cube). Here X 1 is the writer vertex on time step t and 
G* is the colored trinet on time step t. In each of the systems we discuss, we start 
with a cube as our initial network, G° and each structural modification made to 
the network during an update involves replacing a single vertex with a triangle. 
It follows that no vertex will have more than one pair of interlinked neighbors, 
in any network that occurs during the evolution of one of our systems. It follows 
that the type, p(v) of each vertex is well defined. In our systems the current 
type of the writer vertex determines how the writer will move and whether or 
not it will replace its old location with a triangle. 

The following notation will help us to describe motion. For a word w e 
{r, 6, g}* we let v[w] denote the vertex one arrives at by starting at v and 
following the sequence of edges whose colors are described by w. More precisely, 
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for a color c G {r, 6, 17} and vertex u in a colored trinet G let u[c] denote the vertex 
linked to v by an edge of color c. If w is the empty word, we let v[w] — v. For a 
word of the form xc (where x € {r, b, <?}* and c G {r, &, g}) we let v[xc] = v[x] [c]. 



A.l Proof of Theorem Q] 

Let (X 4 , G*) denote the state of the system on time step t, evolving under the 
rules shown in Figure [9] (starting with G° as the cube). Let U(t) denote the 
sequence t/(i)o...?/(t)|[/(t)|_i of vertices in G* that one reaches by starting at 
the writer vertex X 1 and moving along edges with color alternating between red 
and blue. More precisely, Vfc G {0, 1, \U(t) \ — 1} we have 



The length of this sequence, \U(t)\ = 2.min{£; > : X t [(rb) k ] = X 1 }, is the 
number of vertices in the cycle of edges with colors alternating between red and 
blue, that contains the writer. Let 



be the sequence of types of vertices in U(t). In other words, \J(t)\ = \U(t)\ and 
J(t) k = P (U(t) k ) Gt , Vfc g {0, 1, I J(t)\ - 1}. 

Lemma 3 VT > I, we have J(2T) = (bgrO) T+1 and J(2T+ 1) = 00(6yr0) T+1 . 

Proof For T = 1 we have J(2) = bgrObgrO and J(3) = OObgrObgrO, and so the 
result holds in this case. 

Suppose the result holds for some T > 0. 

Now, J(2T+ 1) = 00(6gr0) T+1 , by hypothesis. This means the writer has 
type p(X 2T+1 ) G 2T+i = J(2T+ l)o = on time step 2T + 1. Now, by consulting 
the rules behind the system one can see that updating, from this state, will 
consist of the writer moving along a red edge, and then a blue edge (to the 
position of vertex U(2T + 1)2) and then replacing its previous location with a 
triangle. This implies that the type sequence J(2T + 2), for the next time step, 
can be constructed by taking J(2T + 1) and rotating it two steps to the left 
(to represent to writer's movement) to get (bgrO) T+1 00, and then replacing the 
writers previous position (i.e., vertex corresponding to the penultimate element 



U{t) k 



( 



X t [(rb) k ^ 2 } if k is even, 

X t [r(br)V t ~ 1 V 2 ] if k is odd. 



(2) 



J(t) = J(t) ...J(t)|j(t)|-i = p(U(t) ) G t...p(U(t) m)l _ 1 ) Gt 



(3) 
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in 05rO) T+1 00) with a triangle (bgr) to get J(2T + 2) = (bgrO) T+1 bgrO = 
{bgrO) T + 2 . 

On time step 2T + 2 = 2(T + 1) the writer has type p(X 2(T+1) ) G 2 ( T+i) = 
J(2(T + l))o = b. Now, by consulting the rules behind the system one can see 
that updating, from this state, will consist of the writer moving along a red edge, 
and then a blue edge (to the position of vertex [/(2(T + 1)) 2 ) and then replacing 
its previous location with a triangle. This implies that the type sequence J(2(T+ 
1) + 1), for the next time step, can be constructed by taking J(2(T + 1)) and 
rotating it two steps to the left (to represent to writer's movement) to get 
rO(bgrO) T+1 bg, and then replacing the writers previous position (i.e., vertex 
corresponding to the penultimate element in rO(bgrO) T+1 bg) with a triangle 
(bgr). It follows that the type sequence on the next time step will be J(2(T+1) + 
1) = 00(fegrO) T+1 6<7rO = 00(bgrO) T+2 . The reason for this is that, immediately 
before the update, the writer's previous position (which corresponds to the 
penultimate element in the word rO(bgrO) T+1 bg) is part of a pre-existent triangle 
(together with the vertices corresponding to the last and first elements in the 
word rO(bgrO) T+1 bg), and when the writers previous position is replaced with 
a triangle this pre-existent triangle becomes a square (thus effecting the types 
of the other vertices which used to be part of it). The effect of this is that we 
replace the rightmost b in our word rO(bgrO) T+1 bg with bgr (to represent the 
replacement of the writers previous position with a triangle) and we replace the 
rightmost g and the leftmost r, in our word, with to represent how the types 
of these vertices change when the triangle they used to be part of is expanded 
apart by the addition of the new triangle. These modifications lead to the word 
J(2(T + 1) + 1) = 00{bgrO) T+1 bgrO. 

Now we have shown that our result holds for T = 1, and if our result holds 
for any T > 1 then it holds for T + 1. We can hence prove the result holds 
VT > 1 using induction with T. 

Corollary 4 Vt > 0, the update at time step t consists of the writer moving 
along a red edge and then a blue edge, and replacing its previous position with 
a triangle. 

Proof For t = and t — 1 this result is clearly true. Now suppose t > 2. 

If t is even then 3T > 1 such that t — 2T, and so Lemma [3] implies that 
J{t) = J(2T) = (bgrO) T+1 , and so p(X*) Gt = J(2T) = b, and so (according to 
the rules of our system) the update at time step t consists of the writer moving 
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along a red edge and then a blue edge, and replacing its previous position with 
a triangle. 

If t is odd then 3T > 1 such that t = 2T + 1, and so Lemma [3] implies that 
J(t) = J(2T) = 00(6grO) T+1 , and so p(X*) G t = J(2T) = 0, and so (according 
to the rules of our system) the update at time step t consists of the writer moving 
along a red edge and then a blue edge, and replacing its previous position with 
a triangle. 

We say that n > 2 satisfies property (1) when J(2" — 2) = (bgrO) 2 " \ 
and for each vertex v in G 2 ~ 2 of type p(v)q2™~2 € {r, b} there is a constant 
q(v) £ {0,l,...,2 n -l} such that v = X 2 "- 2 [(rb)^} = U(2 n - 2) 2q{v) (in other 
words, each vertex of type r or b can be reached by repeatedly doing rb hops 
from the writer's current position X 2 ~ 2 ). Here an rb hop consists of moving 
along a red edge, and then a blue edge. 

Suppose that n > 2 satisfies property (1). We will show that this implies 
that n + 1 satisfies property (1), and that Q 2 " +1 ~ 2 can b e obtained by taking 
G 2 ~ 2 and simultaneously replacing each vertex of type r or b with a triangle. 

Now according to Corollary [4] the writer will visit the vertices 

in sequence, over the next 2 n updates, and each such vertex will be replaced with 
a triangle. It follows that the network G 2 "- 2 + 2 ™ — G 2 " +1 ~ 2 present 2™ time 
steps later can be obtained by taking G 2 ~ 2 and simultaneously replacing each 
vertex of type r or b with a triangle. Moreover, each vertex v' of G 2 + ~ 2 which 
is of type r or b will be part of a triangle which was created during an update 
which occurred upon a time step in {(2" - 2), (2™ - 2) + 1, (2" - 2) + 2" - 1}. 
It follows that each such vertex v' will be reachable by doing some number 
q'(v') > of rb hops from the writer's current position X 2 + ~ 2 . Moreover, 
Lemma [5] implies that J(2™ +1 — 2) = (6<?r0) 2 ". This implies that each of our 
vertices v' (with type r or b) can be reached by doing a number of rb hops q'(v') 
such that q'(v) < 2 n+1 — 1. It follows that n + 1 satisfies property (1). 

Now n — 2 clearly satisfies property (1). It follows, by induction with n that 
Vn > 2 we have n satisfies property (1) and G 2 " +1_2 can be obtained by taking 
G 2 ~ 2 and simultaneously replacing each vertex of type r or b with a triangle. 
QED. 
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A. 2 A complete description of the behaviour of the rule 
shown in Figure [9] 

Theorem [I] shows how the system pictured in Figure [9] can be viewed as evolving 
via global replacement operations. Under this system we can describe exactly 
how the network structure and writer position change with time. 

Let N and No denote the sets of positive and non-negative integers respec- 
tively. 

For any pair of non-negative integers x, y £ No let £(x, y) — min{a; + 2y, 2x}. 

We describe the state of a trivalent network automata by listing the vertex 
set, set of red edges, set of blue edges, set of green edges and the position of the 
writer, in that order. 

Suppose T e N : T > 2. Let n(T) = max{m e N : 2 m - 2 < T} and 
t(T) =T- (2"( T ) -2). 

Now let us define the state H{T) = (V{T), E r (T), E b (T), E g (T),W(T)) at 
time T such that: 

1. The set of vertices is V(T) = {-4, -3, .., 2 n{T ^ +1 + 2.t(T) - 1}. 

2. The set of red edges is E r (T) = {{-1, -4}, {-3, -2}}U{{i,i+l mod 2"( T ) +1 + 
2.t(T)} : i G {1, 3, .., 2"( T ) +1 + 2.t{T) - 1}}. 

3. The set of blue edges is E b (T) = {{-4, -3}, {-2, -1}} U {{i,i + 1} : i e 
{0,2,..,2™( T ) +1 +2.i(T)-2}}. 

4. The set of green edges is E g (T) = {{-4, 0}, {-3, C(2™< T ^ 1 , t)}, {-2, ((2 n ( T \ t)}, {-1, ((2 n( - T 1 +1 - 
2"( T )- 1 ,t)}}U{{C(2 Q+2 ./3+2 a ,t),C(2 a+2 ./3+2 a +2 Q+1 ,i)} : a E {0, 1, .., n(T)- 

2}, /3 £ {0, 1, .., 2 n W- a - 1 - 1}} U {{4m + 1,4m + 3} : m € {0, 1, .., t(T) - 
!}}• 

5. The writer position is W(T) = A.t{T) + 1. 

Now we claim, for each T > 2, that H(T) is equivalent to the state of 
S(T) of the system obtained by evolving the rule shown in Figure [9] for T time 
steps, starting from the cube. More precisely, we claim that there is a graph 
isomorphism from the network of H(T) to the network of S(T) that preserves 
the position of the writer and the colors of the edges. Such an isomorphism 
is a bijection from V(T) to the vertex set of the network in state S(T) which 
sends the writer position W(T) to the position of the writer in S(T), and which 
has the property that, a pair of vertices are linked by an edge of a particular 
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color, in H(T), if and only if the images of these vertices under our bijection, 
are linked by an edge of the same color. 

This equivalence of H(T) and S(T) can be seen by induction. It can easily 
be verified in the case where T = 2. Suppose the result holds for some generic 
T > 2 (i.e., suppose that H(T) is equivalent to S(T)). Now we simply have to 
show that the result also holds for T + 1, and we can use induction with T to 
complete the proof. We can show the result holds for T + 1 by showing that 
the state H* obtained by updating H{T) under the system shown in Figure [9] 
is isomorphic to H(T + 1). 

The state H* obtained by updating H(T) is equal to the state obtained 
by taking H(T), and moving the writer along a red edge followed by a blue 
edge, and then performing a structural modification upon the network. This 
structural modification (which yields the network of H* from the network of 
H(T)) consists of the deleting the vertex W{T) (and edgesQ {u bl W{T)} G 
E b {T),{u g ,W{T)} e E g (T),{u r ,W(T)} e E r (T)} that were previously inci- 
dent upon this vertex), and then adding three new vertices, W(T) b , W(T) g , W(T) r , 
together with new blue edges {u b ,W(T) b },{W(T) g ,W(T) r }, and new green 
edges {u g , W(T) g }, {W(T) b , W(T) r }, and new red edges {u r , W(T) r }, {W{T) h , W(T) g }. 
The net effect of these substitutions is to replace W(T) with a triangle upon 
vertices W(T) b , W(T) g , W{T) r . 

We can show H* is equivalent to H(T + 1) by explicitly constructing an 
isomorphism / from H* to H(T + 1) which preserves the position of the writer 
and the colors of the edges. Here / is a bijection from the vertex set of H* to 
the vertex set of H(T + 1) which is defined as follows: 

• For each vertex x of H* such that x < W(T) = 4.t(T) + 1, we have 
f(x) = x. 

• For each vertex x of H* such that x > W(T) — A.t(T) + 1, we have 
f(x)=x + 2. 

• f{W(T) b ) = 4.t(T)+l, f(W(T) g ) = 4t(T)+2 and f(W(T) r ) = 4.t(T)+3. 
A.3 Proof of Theorem [2] 

Let (X*, G*) denote the state of the system on time step t, evolving under the 



rules shown in Figure 12 (starting with G° as the cube). 



Here Uf,,u g and u r denote the vertices which are linked to the writer by edges of colors 
blue, green and red, respectively. 



3G 




B 




Figure 26: A representation of the (essentially one dimensional) dynamics of 
the system illustrated in Figure |12| This picture represents how the relevant 
part of the network changes over successive time steps, like Figure [To} The 
effective rewrite rules are shown at the top right. At the bottom we illustrate 
the symbolic substitutions we shall use to discuss this system algebraically. 



37 



Our method of proof will be to convert the system into a one dimensional 
rewrite system, and then derive a formula for the type p(X t )ct of the writer 
vertex at time t. We can use this information to characterize the number of 
vertices in G*, because the number of vertices in G* only depends on the type the 
writer had on previous time steps (i.e. how much the network grew on previous 
time steps) . The way this system can be converted to a one dimensional rewrite 



system is illustrated in Figure 26 Rows of this figure (reading downwards) show 
the relevant partj^Jof the network on subsequent time steps. 

Let A and B denote the two kinds of vertices illustrated at the bottom of 
Figure 26 We can define these more rigorously, by saying that a vertex v in G* 



of kind A is a vertex which is not part of a triangle, that can be expressed as 
v = X t [r(br) k \^ for some k > and a vertex v of kind B is a vertex which is part 
of a triangle, that can be expressed as v = X t [r(br) k ], for some k > 0. Now we 



can represent our system as a word of symbols that evolves like a tag system 14 
In particular, the word s(t) £ {A, B}* gives the sequence of kinds of vertices 
one visits by starting at X*[r] (the vertex linked to the writer by a red edge) 
and repeatedly doing br hops. More precisely s(t) — s(t)os(t)i...s(t)\ s M\_i and 
Vfc e {0, 1, \s(t)\ - 1} we have 



A i£ p{X*[r{br) k ])a* =0, 
B if p(X*[r(&r) fc ]) G * = g. 



«(*)* = 1 „ ., , vt , ,r (4) 



Here the length of our word is \s(t)\ — minjfc > : X t [(rb) k ] — X*}, which is 
the number of red edges in the alternating rb cycle. 

As the system is updated, the writer modifies the alternating rb cycle and 
moves around it. The way the system gets updated is determined by the type 
of the writer, which is determined by the forms of s(t)\ a (t)\-i and s(t)o (which 
represent the kinds of the vertices to the writer's left and right, respectively). 
In our representation, using s(t), an update corresponds to making s(t + 1) 
by taking s(t), altering the elements at the extreme left or right hand side 
(since these are the parts adjacent to the writer) and also rotating the word (to 
represent the writers movement). 

s(0) = AA 

s(l) = ABA 

s(2) = BAB A 



5 The relevant part of the network G J is the subgraph induced on the cycle with edges of 
colors alternating between red and blue, which contains the writer, with green edges that are 
not part of triangles removed (see Figure 1 10b . 



:-S8 



s(3) = ABAAB 

To see the way s(t) evolves with time, note that Vf > we have: 

1. If s(t) = AxA, for some x € {A,B}* then p(X*) G i = and s(f + 1) = 
xABA. 

2. Ifs(t) = AxB, for some a; G then p{X t ) G t = r and s(m) = BAx. 

3. If s(i) = BxA, for some x € {A, B}* then p(X t ) G t = b and s(i + 1) = 
ABAx. 

These three statements are sufficient to describe the evolution of the system. 
They correspond to the rewrite rules shown at the top right of Figure [26] For 
a word x — XQXi...x m -i let last (a;) = x m —i denote the last element of x. The 
key to our proof is to establish a correspondence between the dynamics of s(t) 
and the dynamics of the binary tag system K(T). Here K(T) denotes the word 
obtained by updating K(0) — 0, a number of times equal to T, under the tag 
system with rules — > 01 and 1 — ^ Oil. More precisely, VT > we have: 

1. If there exists a word x e {0, 1}* such that K(T) = xO then K(T + 1) = 
Olx. 

2. If there exists a word x € {0, 1}* such that K{T) = xl then K(T + 1) = 
Ollx. 

K(0) = 
K{1) = 01 
K(2) = 0110 
K(3) = 01011 

The correspondence between K(T) and s(t) involves the mapping p : {0, 1}* M> 
{A, B}* defined such that for each binary word x £ {0, 1}* we have that p,(x) is 
the word obtained by taking x and replacing each with AB and each 1 with 
A. The correspondence is not trivial. When last(K(T)) = 0, the single step 
K(T) -> K(T + 1) corresponds to two steps, s(t) -> s(t + 1) -> s(t + 2), in the 
s system. When last(if (T)) = 1, the single step K(T) — > K(T + 1) corresponds 
to four steps, s(t) — > s(t + 1) -> s(t + 2) s(t + 3) -» s(t + 4), in the s system. 

VT > let us define the s-time j(T) corresponding with T as 

7(T) = 2 + 2.T + 2.|{t G {0, ...T - 1} : \ast{K{i)) = 1}|. (5) 
Here s(j(T)) represents K(T), as we describe with our next result. 
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Lemma 5 VT > the following equations hold: 



S ( 7 (T)) = Bn{K{T))A (6) 

S ( 7 (T) + 1) = ABA»(K(T)). (7) 
Moreover, if\ast(K{T)) = 1 i/ien </ie following two equations also hold: 

s( 7 (T) + 2) = BAfi(K(T))BA (8) 

s( 7 (T) + 3) = ABAAn{K{T))B. (9) 



Proof For T = we have 7 (T) = 2, #(T) = 0, last(if(T)) = 0, fi(K(T)) = 
AB, s( 7 (T)) = BAB A and s( 7 (T) + 1) = ABAAB, and so our result hold in 
this case. Suppose the result holds for some T > 0. We will show that it holds 
for T+ 1. 

First let us consider the case where last (K{T)) — 0. In this case there must 
exist some x € {0, 1}* such that K (T) — xO, and so we have 

s( 7 (T) + l) - ABA f i(K(T)) (10) 
= ABAp(i)/i(0) (11) 
= ABA[i(x)AB. (12) 

Here, Eq. ( JlQ] ) follows from our assumption that the result holds for T. It 
follows (from the update rules of s) that s( 7 (T) + 2) = BABAfi(x)A. Also, in 
this case, 

7 (T+1) = 2 + 2.(T+l)+2.|{ie{0,..,T}:last(ir(i)) = l}| (13) 
= 7(20 + 2. (14) 

Moreover, K(T + 1) = Qlx, and so we have 

Bn(K{T+l))A = BABA\i{x)A (15) 
= *(7(T + 1)), (16) 

(as we require). It follows (from the update rules of s) that 

s( 7 (T+l) + l) = ABA(i(K{T+l)) (17) 
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(as we require). Moreover, if last(K(T+ 1)) = 1 then there exists a binary word 
y e {0, 1}* such that K(T + 1) = yl. Now 



s( 7 (T+l) + l) = AB^(i/Kl) 
= ABA(i(y)A. 



(18) 
(19) 



Now it follows (from the update rules of s) that 



s( 7 (T+l) + 2) = BA[i(y)ABA 
= BAn(yl)BA 
= BAn(K(T + l))BA, 



(20) 
(21) 
(22) 



(as we require). Now it follows (from the update rules of s) that 



s( 7 (T + 1) + 3) = ABAAfi(K(T + 1))B 



(23) 



(as we require) . So now we have shown that the result holds for T + 1 in this 
case (where last(X(T)) = 0). 

Now let us consider the other case, where last (if (T)) = 1. In this case there 
must exist some x G {0, 1}* such that K(T) = xl. Our assumption that the 
result holds for T gives us s( 7 (T) + 3) = ABAA[i(K(T))B. It follows (from the 
update rules of s) that 



7 (T+1) = 2 + 2.(T+l) + 2.|{ie{0,..,T}:Iast(A'(*)) = l}| (27) 



Hence s( 7 (T+l)) = s( 7 (T)+4) = BABAAfi(x)A. Moreover, K(T+1) = 011a;, 



s( 7 (T) + 4) = BABAA(j,(K(T)) 
= BABAAfj,(x)fj,(l) 
= BABAAfi(x)A. 



(24) 
(25) 
(26) 



In this case, 



BABAAn(x)n(l) 
7 (T) + 4. 



(28) 
(29) 
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and so we have 



Bfi(K(T + l))A = BABAAfi(x)A (30) 

= BABAAfi(x)n(i) (31) 

= s( 7 (T + l)), (32) 

(as we require). Now we can establish that 

s( 1 (T+l) + l)=ABA f i(K(T+l)), (33) 



in exactly the same way we showed Eq. ( 17). In addition, if last(K(T + 1)) = 1 
then there exists a binary word y € {0, 1}* such that K(T + 1) = yl, and we 
can establish that 

s( 7 (T + 1) + 2) = BAn{K(T + l))BA (34) 

and 

s(j(T + l)+3)=ABAA^(K(T + l))B, (35) 



in exactly the same way as we showed Eq. (22) and Eq. (23), respectively. 
So now we have shown that the result holds for T + 1 in this case (where 
last(JC(T)) = 1). So now we have shown that, if the result holds for T then it 
holds for T + 1. Now since the result holds when T = we can use induction 
with T to prove the result holds VT > 0. 

We can use this lemma to form the following corollary, that describes the 
sequence p(X°) G op(X 1 ) G i ... in terms of the forms of K(T) : T > 0. 

Corollary 6 Vt > we have the following: 

1. Ift>0 and t is even then p(X t )ct = b. 

2. Ifte {0, 1} or ift is of the form 

3 + 2.T + 2.|{i € {0,..,T - 1} : last(if(i)) = 1}|, 

for some T > 0, that is such that last (if(T)) = 1 then p{X t )Qt = 0. 

3. Otherwis^we have p(X t )Qt = r. 

6 i.e., when t > 1, t is odd, and there does not exist any T > such that 3 + 2.T + 2.\{i S 
{0,..,T- 1} : last(K(i)) = 1}| 
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Proof Clearly p(X Q ) G o = ppf 1 )^ = 0, so the result holds Vi G {0, 1}. For 
t > 2, let T* be the maximum integer such that j(T*) < t. We must have T* > 
since 7(0) = 2 < t. Also, we must have j(T*) > t — 3, since otherwise (i.e., if 
7(T*) < f - 3) we would have 7(T* + 1) < -y(T*) +4 < f which would contradict 
our assumption that T* is maximum. Now we know j(T*) G {i— 3,i — 2,t — l,t} 
let us consider the different cases. 

If j(T*) — t then t is even and Lemma [5] gives us that s[t) = Bp(K(T*))A and 
so p(X t ) G t = b. 

If 7(T*) = t - 1 then f = 7 (T*) + 1 is odd and Lemma [5] gives us that s(t) = 
s( 7 (T*) + l) = ABAfi(K(T*)) and so last(i^(T*)) = implies p{G t ) X t = r and 
\ast(K(T*)) = 1 implies p(G*)x* = 0. 

If j(T*) = t - 2 then t = 7(T*) + 2 is even and Lemma [5] gives us that s(t) = 
s(Vj(T*) + 2) = BA[i{K(T*))BA and so p(G") xt = 6. 

If 7(T*) = i - 3 then t = j(T*) + 3 is odd and Lemma [5] gives us that s(t) = 
s{ip{T*) + 3) = ABAA(i(K (T*))B and so p(G l ) x t = r. 

These are all the possibilities. So we have V£ > 2 that, p(G t )x t = b if and 
only if i is even. When t > 2 is odd, we must have p(G t )x t G {0, r}. 
Moreover, p{(P)x* = if and only if t = 7 (T*) + 1 where last (If (7 1 *)) = 1. In 
other words, p(G t )x t = if and only if there exists some T > such that t = 
1+t(T) = 3 + 2.T + 2.|{i G {0,..,T-1} : last(if(i)) = 1}| and last (if (T)) = 1. 
Also, if p(G t ) x t i {6,0} then p(G f ) x * = r. 

Now we have established Corollary |6j we can gain a description of the sequence 
p{X°) G op(X r ) G i... by getting a formula for last(if(T)). Our tag system if(T) 
evolves like a slow version of the neighbor independent substitution system Z(n), 
where every element of the word is updated every time step under rules — > 01, 
1 -> 011. 

More precisely, let us define the mapping : {0, 1}* H> {0, 1}* such that 
Vir, G {0, 1}* we have that Q(x) is the word obtained by taking x and replacing 
each with 01 and replacing each 1 with 011. 

Now let Z(Q) = and Vn > let Z(n + 1) = Sl(Z(n)) 

Z(0) = 

Z(l) = 01 

Z{2) = 01011 

Z(3) = 0101101011011. 

For a word, x = XQX\...x m —x let rev(x) = x m ~\x m ^2---Xo denote the reversal 
of the order of characters in the word. 
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Lemma 7 The infinite word last (If (0)) last(if (1))... is equal to the infinite 
word rev(Z(0))rev(Z(l))... 

Proof For a word if (T) = K(T) K(T) 1 ...K(T) lK{T)l _ 1 we have 

K(T + 1) = n(K(T) mT)l _ 1 )K(T) ...K(T) lK{T)l _ 2 (36) 

and 

K(T + 2) = fl(K (T) mT) \-2MK(T) mT)hl )K (T) ...K(T) ]K{T)] _ 3 . (37) 
Following this pattern, one can see that 

last (if (T)) last(if (T + 1))... last (if (T + \K(T)\)) 

is equal to 

K(T) lK{T)l _ 1 K(T) lK{T)l _ 2 ...K(T) = rev(if (T)). 

Moreover, after we have performed |if(T)| updates on if (T), we will get 

if (T + \K(T)\) = ft(K(T) )n(K(T) 1 )...n(K(T) lK{T)l _) = Q(K(T)). (38) 

From this one can see that if (T) evolves like a slow version of the substitution 
system Z(n). When if (T) = Z(n) (as is the case when T = n = 0) we will have 

if (T + \K(T)\) = Sl(K(T)) = n(Z(n)) = Z(n + 1) (39) 

and 

last (if (T)) last (if (T + 1))... last ( if (T + \K(T)\)) = rev(Z(n)). (40) 
Our result can then be proved by induction. 

Lemma 8 Vn > 1 we have Z(n + 1) = Z(n)Z(n)Z(n - 1)...Z(1)1 

Proof Z(2) = 01011 = Z(1)Z(1)1 so the result holds for n = 1. Suppose that 
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the result holds for some n > 1, now 



Z(n + 2) 



Q(Z(n + l)) 

Q(Z(n)Z(n)Z(n - 1)...Z(1)1) 
Z(n + l)Z(n + l)Z(ra)...Z(2)011 
Z(n + l)Z(n + l)Z(n)...Z(2)Z(l)l, 



(41) 
(42) 
(43) 
(44) 



and so the the result holds for n + 1. The result is proved, in general, by 
induction with n. 

For a word x = XQX\...x m ^2%m-\ £ {0,1}* let tog(x) = (1— xo)x\...x m -2(\ — 
x m —l) be the word obtained by altering the characters at x's extreme left and 
right hand sides. 

Let L(x) = x\Xi---x m -\ denote the word obtained by deleting the left most 
element from x. 

Lemma 9 Vm > 1, for any sequence of words a 1 ,a 2 ...,a m £ {01,011} we have 
tog(a 1 )tog(a 2 )...tog(a m ) = J L(a 1 )a 2 ...a m 0. 

Proof This clearly holds when m — 1 since tog(Ol) = 10 = L(01)0 and 
tog(Oll) = 110 = L(011)0. Suppose that the result holds for m > 1. Now 
tog(a 1 )tog(a 2 )...tog(a m )tog(a m+1 ) is equal to L(a 1 )a 2 ...a m tog(a m+1 ), which 
is equal to J L(a 1 )a 2 ...a m 0L(a m+1 )0, which is equal to L(a r )a 2 ...a' n a m+1 as re- 
quired. So the result holds for m + 1. By induction with m the result is proved 
in general. 

Lemma 10 Vn > 1 we have tog(Z (n)) = iev(Z(n)), Z(n)o = and last(Z(n)) = 
1. 

Proof tog(Z(l)) = 10 = (Z(l)), so our result holds when n = 1. Suppose the 
result holds for some n > 1. We write Z(n) = Z(n)o...Z\z(n)\— ii now 



tog(Z(n+l)) 



= tog(fi(Z(n) )...n(^(n)| Z(n) |_i)) 

= tog(01Q(Z(n) 1 )...Q(Z(n)| X(n)h2 )011) 

= lin(Z(n) 1 )...n(Z(n)| Z(n) |_ 2 )010, 



(45) 
(46) 
(47) 
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also, 



rev(Z(n + l)) = rev(n(Z(n)| Z(n) |_ 1 ))...rev(n(Z(n) )) (48) 

= 110rev(n(Z(n)| Z(n) |_ 3 ))...rev(n(Z(n)i))10 (49) 

= 110rev(n(Z(n) 1 ))...rev(n(Z(n)| Z(n) |_ 2 ))10 (50) 

= 110tog(O(Z(n) 1 ))...tog(O(Z(n)| Z( „ ) |_ 2 ))10 (51) 

= 110i(n(Z(n)i))n(Z(n) 2 )...n(Z(n)| Z(n) |_ 2 )010 (52) 

= lin(Z(n)i)n(Z(n) 2 ...n(Z(n)|z(„)|-a)010 (53) 



The reason Eq. (49) implies Eq. (501 is that rcv(Z(rt)) = tog(Z(n)) implies 

Z{n)x...Z(ri)\ Z {n)\-2 = Z(n) ]z{n)] _ 2 ...Z{n)i 
is palindromic, which implies 

rev(0(Z(n)| Z(n) |_ 2 ))... rev(n(Z(n) 1 )) = rev(0(Z(n) 1 ))...re V (^Z(n)| Z(n) |_ 2 )). 



The reason Eq. (501 implies Eq. |5T| ) is that we always have f2(Z(n)j) € 
{01,011}, and so rev(fi(Z(n)j)) = tog(f2(Z(n)j)). We can use Lemma [9] to 



derive Eq. (52) from Eq. (51 ), and Eq. (|53j follows because the left most entry 



in Q,{Z{n)x) is 0. 



Now, since Eq. (53) is identical to Eq. (471, we have proved that our result 



holds for ft + 1. The fact that the result holds Vn > 1 then follows by induction. 

Now we can state our critical result. Recall that last(i£"(0)) last(if (1))... is 
the infinite sequence consisting of the last elements of the words K(0), K(l)... 
(evolving under our tag system). 

Lemma 11 last(if(0)) last (if (1))... = ]im n ^. 00 Z(n) 

Proof Lemma [8] gives us that Vn > 0, 

Z{n + 1) = Z{n)Z{n)Z{n - 1)...Z(1)1, (54) 

and so 



rev(Z(n + 1)) = 1 rev(Z(l)) rev(Z(2))... rev(Z(n)) rev(Z(n)), (55) 
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Now, Lemma 
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implies rev(Z(k)) = tog(Z(k)), Vfc > 0. So now we have 



Z(n + 1) = rev(rev(Z(n + 1))) (56) 

= tog(rev(Z(n+l))) (57) 

= 0rcv(Z(l))rev(Z(2))...rev(Z(n))i?(rev(Z(n)))l, (58) 

where R(x) denotes the word obtained by deleting the right most clement from 
x. Now, taking the limit as n — > oo we have 

lim (Z(n)) = lim(Z(7i+l)) (59) 

n— >oo n—too 

= 0rev(Z(l))rev(Z(2))... (60) 

= rcv(Z(0))rev(Z(l))rev(Z(2))... (61) 

= last(2f(0))last(#(l))last(isr(2))... (62) 

Here Lemma [7] gives us that Eq. (61) implies Eq. (62). 



The above result (Lemma 11 ) is important because a formula is known 20 
for the infinite sequence of lim n _ >00 (Z(n)) (which we have just shown to be 
equal to Iast(if(Q)) last(i^(l))...). Using this formula we have Vi > 0, that 

last(K (i)) = ( lim (Z(ri)j) = Hi + 1U\ - [i.<t>\ - 1, (63) 

where <f> = is the golden ratio. Other results are also known for this 

sequence. For T > the index of the Tth 1 in the sequence lim„_ ) . 00 (Z(n)) is 



well known 21 to be given by the formula [(T + l)<f>\ ■ In other words, 
{ T > : ( lim (Z(n))) = 1 1 = {T > : last(K(T)) = 1} = { [{n+l)(j)\ :n>0} 

[ Vri— >oo / T J 

(64) 

Suppose h > 2. Now, accor ding to Corollary § p(X h ) Qh = if and only if 
3T > such that h = 3 + 2T + 2|{i G {0, 1, T - 1} : last (iT(i)) = 1}| and 
last(fc(T)) = 1. 

Now since T > 0, last(fc(T)) = 1 if and only if 3n > : + 1)0J = T we 
have the following: 



47 



p(X h ) Gh = if and only if 3n > such that h = 3 + 2|_(n + 1)^J + 2.9, where: 



9 = |{ze{O,l,...,L(n+l)0J-l}:last(if (*)) = !}! 



(65) 



|({0, 1, L(n + 1)0J - 1} n {L(o + l)^j , L(l + 1)0J , (66) 



= |{L(O + l)0j,L(l + l)0j,...,L((n-l) + l^j}| 
= n. 

Hence we have 

{h > 2 : p(X h ) G H = 0} = {3 + 2[(n + 1)0J + 2.n : n > 0}. 

Now, the number of vertices |G*|, in G* will be equal to 

|G*| - 8+2|{z G {0, ...,i-l} : - 6}|+2|{i G {0, .„,<-!} : p{X% 



(67) 
(68) 



(69) 



= 0}| 
(70) 



because the number of vertices only changes when the writer is of type or 6, 
and it increases by 2 each such time. Now Corollary [6] implies 



\{i€{0,...,t-l}:p(X% i =b}\ = 



t+1 



-1 + 4, 



(71) 



(where Si j is the Kronecker Delta). 
Also, Vt > 5, we have: 

|{t G {0, ...,t- 1} : p(X%t = 0}| = 2 + |{t G {2,3..., t - 1} : p{X i ) & = 0}| 

(72) 

\{i G {0,...,t-l} : p{X l ) G , = 0}| = 2 + \{n > : 3 + 2L(n + l)0j + 2.n < *}l> 

(73) 



Where we use use Eq. ( 69 ) to infer Eq. ( 73 ) from Eq. ( 72 ). Now Vn > Vt > 5 



the following four inequalities are equivalent: 

'S + 2[(n + l)(f>\ +2.n<t 



l(n + l)4>\ < 



-(n + 1) 



(74) 
(75) 
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{n + l)(j) 2 < 



n < 





t 




2 



(76) 



(77) 



The equivalence of inequalities ( 74 ) and ( 75 1 can be derived using the fact that 
2a < b a < [^pj , for integers a and 6. The equivalence of inequalities (751 
and ( 76 1 can be derived using the fact that [r\ < m r < m for a real r and 
an integer m, and the fact that 1 + cf> — (fx 2 . 



It follows that when t > 5 we can rewrite Eq. (|73|) as 
\{iE{0,...,t-l}:p(X l ) G > =0}| 



= 2 
= 2 

= 2 



n > : n < 
0,1. 



1 t 





t 


1 " 














2 


02 





1 



1 



8) 
(79) 
(80) 



Now by substituting Eq. (711 and Eq. (80 1 into Eq. (70), we get that when 
t > 5, the number of vertices in G t is 

t + 1 



|G £ 



8 



1 



2 2 



(81) 



as required. It can also be checked that Eq. (81 1 gives the correct answer for 
each t e {0,1,2,3,4,5}. Q.E.D. 
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